Optics of a Faraday-active Mie Sphere 

D. Lacoste, B. A. Van Tiggelen 
Laboratoire de Physique et Modelisation des Milieux Condenses, Maison des Magisteres, 

B.P. 166, 38042 Grenoble Cedex 09, France 
G.L.J. A Rikken, and A. Sparenberg 
Grenoble High Magnetic Field Laboratory, Max-Planck Institut 
fiir Festkorperforschung/C.N.R.S, B.P. 166, 38042 Grenoble Cedex 9, France 



Abstract 



We present an exact calculation for the scattering of light from a single sphere 
made of a Faraday-active material, into first order of the external magnetic field. 
When the size of the sphere is small compared to the wavelength the known T-matrix 
for a magneto- active Rayleigh scatterer is found. We address the issue whether or not 
there is a so called Photonic Hall Effect - a magneto-transverse anisotropy in light 
scattering - for one Mie scatterer. In the limit of geometrical optics, we compare our 
results to the Faraday effect in a Fabry-Perot etalon. 

PACS: 42.25.Fx, 78.20.Ls 
1. Introduction 

Several reasons exist why one wishes to understand light scattering from a dielectric sphere 
made of magneto-active material. Single scattering is the building block for multiple scat- 
tering. Many experiments have been done with diffuse light in a magnetic field&i. Though 
qualitatively very useful, it turns out that a theory for point-like scatterers in a magnetic 
field, as first developed by MacKintosh and JohnB and later refined by Van Tiggelen, May- 



nard and Nieuwenhuizem, does not always describe observations quantitatively, for the 




evident reason that experiments do not contain "small" scatterers. This paper addresses 
light scattering from one sphere of any size in a homogeneous magnetic field. 

The model of Rayleigh scatterers has been successfully used to describe specific proper- 
ties of multiple light scattering in magnetic fields such as coherent backscattering and the 
Photonic Hall Effect . In the following section we present the perturbation approach on 
which our work is based, which will allow us in section 3 to compute the T-matrix for this 
problem. With this tool, we will be able to answer in section 4 the issue concerning the 
Photonic Hall Effect . 

2. Perturbation theory 

We consider the light scattered by one dielectric sphere made of a Faraday- active medium 
embedded in an isotropic medium with no magneto-optical properties, using a perturbative 
approach to first order in magnetic field. 

In this paper we set cq = 1. In a magnetic field, the refractive index of the sphere is a 
tensor of rank two. It depends on the distance to the center of the sphere r, which is given 
a radius a via the Heaviside function 0(r — a), that equals to 1 inside the sphere and 
outside, 



In this expression, Eq = is the value of the normal isotropic dielectric constant of the 
sphere of relative index of refraction m ( which can be complex ) and Sp = 2mVo/c(j is a 
coupling parameter associated with the amplitude of the Faraday rotation ( Vq being the 
Verdet constant and uj the frequency ). We introduce the antisymmetric hermitian tensor 
$jj = ieijkBk- Except for ef, the Mie solution depends on the dimensionless size parameters 
X = ka with k = uj/cq and y = mx. In this paper, we restrict ourselves to non-absorbing 
media, so that m and ep are real- valued. 
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Noting that the Helmholtz equation is formally analogous to a Schrodinger equation with 
potential V(r,ti;) = [I — £:(B,r)] liP' and energy o;^, the T-matrix is given by the following 
Born serieJ^: 

T(B, r, uj) = V(r, u) + V(r, u) ■ Go ■ V(r, u) + ... (2) 

Here Go{uJ,p) = l/{uj^l — p^Ap) is the free Helmholtz Green's function and {Ap)ij = 
Sij — PiPj/p'^- We will call T° the part of T that is independent of the magnetic field and 
the part of the T-matrix linear in B. It follows from Eq. (§) that: 

= (l-£o)£peTO-$- (I + Go-T°). (3) 

We need to introduce the unperturbed eigenfunctions \E'^j(. (r) of the conventional Mie prob- 
lem. These eigenfunctions represent the electric field at the point r for an incident plane 
wave I (T, k > along the direction k with an helicity a. This eigenfunction is "outgoing" 
for \t'^,^ and "ingoing" for \E'~j^ according to the definition of the outgoing and ingoing free 
Helmholtz Green's function: 

|vl/±^ >= (I + G±-T°)|a,k> . 

For our free Helmholtz Green's function, this implieJ^: 

^;;-k(r) = (-l)^+'^^l.,k(r). (4) 

We denote by k the incident direction and k' the scattered direction. With these nota- 
tions, it is possible to obtain from Eq.(^, 

TL,k'.' =eFu;'< I e $ I > . (5) 

This equation can also be obtained from standard first-order Rayleigh-Schrodinger pertur- 
bation theoryi. 

Two important symmetry relations are found to be satisfied by our T-matrix. The first 
one is the parity and the second one is the reciprocity: 



T-k-(j,-k'-cr'(B) — Tko-,k'o-'(B) 



(6) 



-k'cr'.-ko-l 



-B) = Tk.,k'a'(B) (7) 



3. T matrix for Mie scattering 

In order to separate the radial and the angular contributions in Eq. we expand the Mie 
eigenfunction \E'^j^ in the basis of the vector spherical harmonicsi, 

(r) = - ^-^^^Y^;, (f ) //,,(r) Yf^ik) ■ x'„. (8) 

In this definition, p = kmr, the x'a the eigenvectors of the spin operator in the circular 
basis associated with the direction k, and implicit summation over the repeated indices 
J, M, A and A' has been assumed ( the indices A and A' for the components of the field 
can take three values, one being longitudinal and two of them being perpendicular to the 
direction of propagation). 

Because of the presence of the function in Eq. @, we only have to consider the field 
inside the sphere, whose main features are contained in the radial function fx>\{f^) ■ This 
matrix fy-^ is known in terms of the transmission coefficients cj and dj of Ref.i, and of the 
Ricatti-Bessel function uj{p). We found the following expression: 



^ -tu'j{p)cj ^ 



f"'(p) = m 



-tu'jip)cj^J{J+l)/p 
y Uj{p)dj J 

The three vectors Xa - similar to the x'a but associated with the 2-axis - are a convenient 
basis for this problem since they are the eigenvectors of the operator $ with eigenvalue 
— cr, provided we choose the z-axis along B, which we will do in what follows. Eq. (|^) is 
simplified by this choice and the angular integration leads eventually to: 



where Q is the matrix 



Q(J,M) = -MB 



1 



^/J(J+l) 














(9) 



j(j+i) / 

and B being the absolute value of the applied magnetic field. The linear dependence on the 
magnetic quantum number M can be expected for an effect like the Faraday rotation, affect- 
ing left and right circularly polarized light in an opposite way similar to Zeeman splitting. 
The radial integration can be done using a method developed by Bott et al.i. 

IGtt 



rpl 

-■-Ick' 



J,M 



wj:i-M) Cj Y 1,,^ (&) Y l;,^ + vj y|,^ m yJ:^ 



(10) 



with the dimensionless parameter: 



W = VoBX, 



and the coefficients: 



C 



J 



'A 



J 



V 



J 



A*j Aj 



(11) 



(12) 



with Aj{%j) = u'j{y)/uj{y) and B the amplitude of the magnetic field directed along the 
unit vector B . Absorption in the sphere is still allowed. We will consider the limiting case 
of a perfect dielectric sphere with no absorption ( '^m{m) ^ ). Using I'Hospital's rule in 
Eqs. (|lT]) and (p!2|), we obtain immediately for this case: 

Aj J {J + oo) 



J 



oo^ A^ 



V 



J 



( A J J{J + oo) 



J{J + oo) 



t 



oo + A% 



J 



(13) 



(14) 
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A. T-matrix without magnetic field 

For future use, we need the on-shell T-matrix of the conventional Mie-problemi. It is given 
by a formula analogous to Eq. (p!OD where Cj and T>j are replaced by the Mie coefficients 
aj and 6j, and with M = 1. Because of rotational invariance of the scatterer, it's clear that 
the final result only depends on the scattering angle Q which is the angle between k and k' ( 
see Fig. |1] ). Therefore, we get in the circular basis ( associated with the indices o and a' ): 

= ^ E jjJTl) ^""^ + ""^'^^^ + ^^'rjAcosO)] . (15) 

In this formula, the polynomials ttj^m and tjm are defined in terms of the Legendre poly- 
nomials Pj ^ byi 

7rj,M(cose) = (cos^) , tj^m{cos9) = ^Pf{cos9). (16) 

sm u do 



B. Tj^,^, when k / k' 

It remains to express the vector spherical harmonics in Eq. (p!0| ) in terms of the natural angles 
of the problem in the presence of a magnetic field. The latter breaks rotational invariance. 
Because is linear in B, it can be constructed by considering three special cases for the 
direction of B. If k 7^ k', we can decompose the unit vector B in the non-orthogonal but 
complete basis of k, k' and g = k x k'/|k x k'|, and this results in: 

, (B-k)(k-k') -B-k' -, (B-k')(k-k') -B-k , , . ^ 
Tkk' = ^ ^ + ^ t + B ■ g 17 

(k-k')2-l ^='^ (k-kO^-l ^='' ^ ' ^ ' 

with 

Tl,{h = k) = — 5] ll±±-(-a) {Cj + aa'Vj) [vr^.oolcos^) + a<y'rj,^{cos9)] , (18) 



j> 



\J{J + 1] 



Tj,,(B = k') = — E 4tTTT(-^') + ^^'^-7) [^^-oo(cos^) + aaV^,oo(cos( 



(19) 



UJ j^, J(J+1) ^ '{J + M)V ^ ^' 



J>M>0 



(20) 



6 



C. Particular case for T for forward scattering 

The treatment in (3.2) becomes degenerate when k and k' are parallel, i.e. in forward 
scattering. In this case, B can still be expressed in a basis made of k and of two vectors 
perpendicular to k. The contribution of these last two vectors must have the same form as 
in Eq. (^Dj) for = 0. Hence there is no such contribution and we find 

TLm' = (B • k)(-a) — E(2^ + 1) (Cj + Vj). (21) 

^ j>i 

In Fig. we plotted real and imaginary part of this expression in units of as a function 
of the size parameter x for a = — 1 and B = k. The forward-scattering amplitude has an 
important application in inhomogeneous media, namely as the complex average dielectric 
constant. 



D. Optical Theorem 



We will check our formula on energy conservation as expressed by the Optical TheoremE 



u yJ (4"-) 

To first order in magnetic field, the r.h.s of this equation equals 



If we assume that B || k, we can compute this using Eqs. (|18|) , (|15|), and the following 



orthogonality relations for the polynomials ttji and rji ( which we denote as ttj and tj )a : 

d{cos9) [tt J {cos 9)tk {cos 9) + tj{cos9)ttk{cos9)] = 0, 

r 2j2fJ + l)2 

j d{cos9) [ttj{cos9)ttk{cos9) + tj {cos 9)tk {cos 9)] = + I ^^^^ 

The l.h.s of Eq. (P^D, is obtained from Eq. (pi]). The Optical Theorem provides us a 
relation between Mie coefficients, which we can actually prove from their definitions: 
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Ke(a}4^) = 53m(4), (24) 
and 

^e{b*j(fj'^) = Qm{d]). (25) 



4. Magneto-transverse Scattering 

From the knowledge of the matrix T^, we can compute how the magnetic field affects the 
differential scattering cross section summed over polarizations as a function of the scattering 
angle. Only the diagonal part of this matrix in a basis of linear polarization will affect 
the scattering cross section since we consider only terms to first order in magnetic field. 
Therefore only the third contribution of Eq. ( p!7D plays a role, which means that the effect 
will be maximum when B || k x k' ( the typical Hall geometry ). Indeed symmetry implies 
that the magneto-scattering cross section should be the product of det(B, k, k') and of some 
function that entirely depends on the scattering angle 6. 

We choose to normalize this magneto scattering cross section by the total scattering cross 
section in the absence of magnetic field, 



- aa' 



- sin (f)F{e) (26) 



C J,^d cos 9 E.a' IT^P 

The so called Photonic Hall effect (P.H.E) is a manifestation of a magnetically induced 
transverse current in the light transport which has similarities to the Hall effect known for 
the transport of electrons. In an experiment on the P.H.E, one measures the difference in 
scattered light from two opposite directions both perpendicular to the incident direction of 
light and to the applied magnetic fieldi. The P.H.E is a manifestation of the anisotropy of 
light scattering due to a magnetic field in the regime of multiple light scattering. 

Although experiments dealt with multiple scattering so far, it is interesting to see if a 
net magneto-transverse scattering persists for only one single scatterer. In Figs. the 
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magnetic field is perpendicular to the plane of the figure, the incident light is along the x- 
axis. A typical measurement of the magneto-transverse scattered light is therefore associated 
with the projection of the curve onto the y-axis, which we define as the magneto-transverse 
direction. 

5. Magneto-transverse scattering as a function of the size parameter 

Quantitatively, the transverse light current difference is associated with a summation of 
the magneto scattering cross section over outgoing wavevectors and normalized to the total 
transverse light current. A schematic view of the geometry is displayed on Fig. |l|. In our 
notation this is 

^ Iuv{B) - I,oUB) ^ 2J^dcPJ^dcos9 singsin0 E..' ^ejT^^^T^:,) 
I^p{B = 0) + I,on,n{B = 0) d(P J^d cos 9 sm 9 sin (PEa^'lTW ' 

The factor sin 6' sin represents a projection onto the magneto transverse direction B x k 
which is necessary since we are interested in the magneto-transverse light flux. In Fig. |^ 
we plotted this contribution as a function of the size parameter x for an index of refraction 
m = 1.0946 ( the value in Ref.i ). Note the change of sign beyond x = 1.7, for which we do 
not have any simple explanation so far. In the range of small size parameter, t] exhibits an 
x^ dependence. 

A. Rayleigh scatterers 

For Rayleigh scatterers, the formulas (|1^ to (^TJ) simplify dramatically because one only 
needs to consider the first partial wave of J = 1 and the first terms in a development in 
powers of y ( since y < 1 ). From Eqs. and d^), we find : Coo = -Gt^/$^(G + i^^^)^ 
and = — t^/Ay^^ so that we can keep only Coo and drop "Doo as a first approximation. 
Adding all contributions of Eqs. ( [T7| ) and (p!5D, and changing from circular basis to linear 
basis of polarization we find 
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k.k' 



( tok • k' + itiB ■ (k X k') itiB ■ k ^ 
y —itiG ■ k' to J 



(28) 



with to = —Qii^aX/oj the conventional Rayleigh T-matrix and ti = GCooW/cij. This form 
agrees with the Rayleigh point-like scatterer model discussed in Ref.i. We note that Eqs. 
(|I^) and (|TP|) give off-diagonal contributions in Eq. (^) whereas Eq. (^) gives a diagonal 
contribution. This is a general feature that persists also beyond the regime of Rayleigh 
scatterers. 

For a Rayleigh scatterer the magneto cross section of Fig. ^ exhibits symmetry, since 
the positive and negative lobes of the curve are of the same size but of opposite sign. Hence 
no net magneto-transverse scattering exists for one Rayleigh scatterer. In fact Eq. (pH]) 
provides the following expression for F{6): 



'Rayleigh(^) 47^2(^2 ^ 2)2 



cos e sin e. (29) 



As the size of the sphere gets bigger, the magneto corrections become asymmetric, as seen 
in Figure ^. When the size is further increased, new lobes start to appear in the magnetic 
cross section corresponding to higher spherical harmonics. These lobes do seem to have a 
net magneto-transverse scattering. 

One single Rayleigh scatterer does not induce a net magneto-transverse flux. It is in- 
structive to consider the next simplest case, namely two Rayleigh scatterers positioned at 
ri and r2. If their mutual separation well exceeds the wavelength, it is easy to show that 
the collective cross-section simply equals the one-particle cross-section multiplied by an 
interference factor ^(k, k') = |exp(i(k — k') ■ ri) + exp(i(k — k') ■ r2)f. This interference 
factor changes the angular profile of the scattering cross-section and makes sure that a net 
magneto-transverse flux remains. The estimate for two Rayleigh particles with an incident 
wave vector along the inter-particle axis, is found to be: 

VoB 3 /sin(A;ri2)y . , . 

v-^^ —T- — if > 1 (30) 



k \ kr 



12 
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This simple model suggests that the "Photonic Hall Effect" is in fact a phenomenon gener- 
ated by interference of different light paths. In Fig. ^ we show how differential cross-sections 
of two particles change in a magnetic field. More scattering is now directed into the forward 
direction and as a result the cancellation of the net magneto-transverse flux in Fig. |^ is 
removed. One Mie sphere mimics qualitatively this simple model and should on the basis of 
the principle outlined above exhibit a magneto-transverse current. The model also suggests 
that the regime of Rayleigh-Gans scattering! - the Born approximation for one sphere but 
contrary to Rayleigh scattering still allowing interferences of different scattering events - 
should exhibit a net magneto-transverse flux. Indeed, explicit calculations in this regime 
confirm this statement, with ~ x^, independent of the index of refraction m of the sphere. 



B. Geometrical Optics 

In the regime of large size parameters, the Mie solution can be obtained from Ray Optics. 
Apart from Fraunhofer diffraction ( persisting for any finite geometry ), the Mie solution 
for a ray with central impact should be equivalent to the one for a slab geometry. In our 
magnetic-optic approach, this means to study Faraday rotation in a Fabry- Perot cavity. This 
model is of special interest because the Fabry- Perot cavity is known to enhance the Faraday 
rotation0'0, i.e the rotation is additive in the total traversed path length as opposed to the 
case of rotary power. 

A ray with central impact is characterized by J = 1 in Mie theory. We assume that 
X ^ 1 and y ^ 1 which allows some simplification in the expression of the Mie coefficients. 
We find for the Mie coefficients Ci and di the following behavior: 

~ (m + 1) (1 + re-^'y) ' 

^1 = (m + 1) (1 - re-^'y) ' ^^^^ 
In this formula, r = (m — l)/(m + l) is the complex Fresnel reflection coefficient. Putting 
this expression into Eqs. (|T3|) , ([l^ ) and (^T]) , we can compute the exact behavior of the T- 
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matrix in the forward direction. We note T^catt P^^^ of T'' due to scattering only, which 
is obtained from Eq. (pr5| ) by replacing the Mie coefficients aj and bj by aj — ^ and 6j — |, 
since the terms | are associated with the Fraunhoffer diffraction, which does not exist for 
the slab geometry. 

Since for a small perturbation, we have 

T = TL. + T^-TL.e^VTL.., (32) 

we see that the change in phase 5(j) due to the magnetic field is in fact related to the imaginary 
part of T^/T°^^j^. This change in phase can be interpreted as the Faraday effect. 
From Eq. (|3l|) , we find in the basis of circular polarization: 
/ \ 

' scatt , 



Qmij;^] =6<Pi-a)B-k6^^>, (33) 



with 

66 = 2aVoB^—- (34) 

^ " l-i?[l + A^sin(et)e] ^ ' 

with Ai = A7^/(oo — TZ)'^ and the refiectivity R = r^. We note that the quantity (— cr)B • k 

is conserved for a given ray, which generates the accumulation of the Faraday rotation. The 

function 6(j) tends to 2aVo-B as i? — > 0, since it represents the normal Faraday rotation in an 

isotropic medium of length 2a, as it should be for our geometry. When R is large, two new 

factors come into play: {1 + R)/{1 — R) which is the maximum gain factor of the Faraday 

rotation due the multiple interference in the Fabry-Perot cavity and 

Am = ^ 

cx) + A<sin(Gt)e 

which is an Airy function of width 4/ y/jA. The finesse of the cavity is then = ny/jA/ G. 
At resonance, the Faraday rotation is maximally amplified - assuming no losses - relative 
to single-path Faraday rotation0'0. We stress that one needs 50 <^ 1, in order for Eq. 
(p^) to apply. The Faraday rotation has the effect of splitting each transmission peak in 
the Fabry-Perot cavity into two peaks of smaller amplitude both associated with a different 
state of helicity. 
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The amplification of tlie Faraday rotation is a consequence of the amphfied path length 
of the light. In other words, the Faraday rotation measures the time of interaction of the 
light with the magnetic field0. This time is found to be the dwell time r of the light in the 
cavity for this ID problem . The change in phase follows the simple relation: 

5<j) = VoB^ (35) 
m 

where co/m is recognized as the speed of light in the sphere. 

The dwell time of the light in the cavity varies between a maximum value of 

rZel = (1 + m')a/co, 

and a minimum value of 

rSz = 4mV(l + m2)a/co. 

These tj^ical oscillations are visible in the plot of the change of phase 5(j) of Fig. ^ 



6. Summary and Outlook 

In this paper we addressed the Faraday effect inside a dielectric sphere. We have shown that 
this theory is consistent with former results concerning the predictions of the light scattered 
by Rayleigh scatterers in a magnetic field. It is possible to get from this perturbative theory 
quantitative predictions concerning the Photonic Hall Effect for one single Mie sphere, such 
as the scattering cross section, the dependence on the size parameter or on the index of 
refraction. 

We will start experiments addressing single Mie scattering in a magnetic field. A second 
challenge is to implement our Mie solution into a multiple scattering theory. 
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FIGURES 

Fig. 1. Schematic view of the magneto-scattering geometry. Generally, denotes the angle 
between incident and outgoing wave vectors; <j) is tlic azymuthal angle in the plane of the magnetic 
field and the y-axis. The latter is by construction the magneto-transverse direction defined as the 
direction perpendicular to both magnetic field and incident wave vector. The angle a coincides 
with the angle 6 in the special but relevant case that the incident vector is normal to the magnetic 
field. 

Fig. 2. The solid line is the real part and the dashed line the imaginary part of the magneto 
forward scattering matrix Tjj.j^ in the circular basis of polarization, plotted versus the size parameter 
X for an index of refraction of m = 1.33 in units oiW = VqBX. 

Fig. 3. Magneto-transverse scattering cross section F{6) for a Rayleigh scatterer with index of 
refraction m = 1.1 and size parameter x = 0.1. The solid line is a positive correction and the 
points denote a negative correction. The curve has been normalized by the parameter W . No net 
magneto-transverse scattering is expected in this case because the projection onto the y-axis of 

these corrections cancel. 

Fig. 4. Magneto-transverse scattering cross section F{9) for a Mie scatterer of size parameter 
x = 5 and of index of refraction m = 1.1. The curve has been normalized by the parameter W . 
The solid line is for positive correction and the points denote a negative correction. In this case 
a net magneto-transverse scattering is expected because the projection onto the y-axis of these 
corrections do not cancel. 

Fig. 5. Normalized magneto-transverse light current 77 as a function of the size parameter x for 
an index of refraction of m = 1.0946. The curve is displayed in units of W . 

Fig. 6. Magneto-cross section for two Rayleigh scatterers each of size parameter fca = 0.1 and 
separated by a distance corresponding to a size parameter of kr\2 = 5. In this case the enhanced 
forward scattering leads also to a net magneto-transverse current along the vertical axis. 
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Fig. 7. Magnetically induced change of phase d4> - similar to Fabry-Perot modes of a cavity - 
a function of size parameter x for the partial wave of J = 1 - "central impact" -. The curve is 
m = 10, and has been normalized by the value 2aVoB. 
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